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Abstract 

We generalize the concept of ascending runs from permutations to 
Cayley trees and mappings. A combinatorial decomposition of the cor¬ 
responding functional digraph allows us to obtain exact enumeration for¬ 
mulas that show an interesting connection to the Stirling numbers of the 
second kind. Using a bijective proof of Cayley’s formula, we can show why 
the numbers of trees with a given number of ascending runs are directly 
linked to the corresponding numbers for mappings. Moreover, analytic 
tools allow a characterization of the limiting distribution of the random 
variable counting the number of ascending runs in a random mapping. 
As for permutations, this is a normal distribution with linear mean and 
variance. 

Keywords: random mappings; Cayley trees; ascending runs; exact enu¬ 
meration; limiting distribution; Injections 


1 Introduction 

Random n-mappings, or simply random functions from the set of integers [n] := 
{1,2,..., tt.} into itself, appear in various applications: for example in the birth¬ 
day paradox, the coupon collector problem and occupancy problems, just to 
name a few. Structural properties of the functional digraphs of random map¬ 
pings have widely been studied, see e.g. the work of Arney and Bender [I], 
Kolchin 114], Flajolet and Odlyzko |9j. For instance, it is well known that the 
expected number of connected components in a random n-mapping is asymp¬ 
totically 1/2 • log(n), the expected number of cyclic nodes is yJ-Kn/2 and the 
expected number of terminal nodes, i.e., nodes with no preimages, is e~ l n. 

In functional digraphs corresponding to random mappings the nodes’ labels 
play an important role and thus it is somewhat surprising that occurrences of 
label patterns have, so far, not received much attention. The simplest of all 
label patterns in mappings are ascents, i.e., nodes with label i for which it holds 

‘Both authors were supported by the Austrian Science Foundation FWF, grant P25337- 
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that i < f{i) where f(i) denotes the image of i. In a random n-mapping, it 
is clear that the probability of i being an ascent is simply equal to {n — i)/n. 
Ascents in random mappings are thus rather easily dealt with. Egecioglu and 
Remmel [7| showed how results on ascents in mappings can be translated to 
results for ascents in Cayley trees by providing a family of weight preserving 
bijections. As shown by Clark [3], these results can be used to provide central 
limit theorems for the number of ascents in random mappings and in random 
Cayley trees. In this context we also mention the recent paper of Seo and Shin 
[18] , where the size of the maximal subtree without ascents containing the root 
node has been studied for Cayley trees. Some recent research concerned with 
patterns formed by the labels in a functional digraph has been performed in 
ubi. where alternating mappings have been studied. These are a generalization 
of the concept of alternating permutations to mappings. They can be defined as 
those mappings for which every iteration orbit forms an alternating sequence. 
Alternating mappings can thus also be seen as mappings that do not contain 
two consecutive ascents or descents. Also, we would like to mention the PhD 
thesis of Okoth m who has very recently studied local extrema in trees (called 
sources and sinks there). His studies also led to results for the corresponding 
quantities in mappings. 

In this paper we perform an analysis of a fundamental label pattern in random 
mappings. Namely, we generalize the notion of ascending runs from permuta¬ 
tions to mappings. In the corresponding functional digraphs, ascending runs are 
maximal ascending paths. Our goal is to enumerate mappings by their size and 
by their number of ascending runs and to characterize the typical behaviour of 
this label pattern in random mappings. 

This paper is organized as follows: We start by briefly recalling the basic 
definitions that will be used. In Section[3j we present a bijective proof of Cayley’s 
formula that will turn out to be useful for the study of ascending runs in Cayley 
trees and mappings. Section[4]is devoted to the main topic of this paper, namely 
to the exact and asymptotic study of ascending runs in mappings. This leads to 
results for Cayley trees as well. We use a combinatorial decomposition of Cayley 
trees and mappings according to their smallest node. This leads to a recurrence 
relation for the studied numbers that can be translated into a partial differential 
equation for the corresponding generating functions. Solving this PDE allows for 
extraction of coefficients and for a characterization of the limiting distribution 
of the number of ascending runs in a random mapping or tree. The bijection 
presented in Section [3] will allow to explain a simple connection between the 
results for trees and for mappings. Finally, we close this paper by setting our 
results in relation with the corresponding results that have been obtained for 
permutations and pointing out directions for future research. 


2 Preliminaries 

The combinatorial objects of interest in this paper are Cayley trees and map¬ 
pings. These are defined as follows: 

Definition 1. A Cayley tree of size n is a unordered rooted labelled tree with n 
nodes. In a labelled tree of size n every node carries a distinct integer from the 
set [n] as a label. 
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Figure 1: Functional digraph of a 19-mapping. 


We always consider the edges in a Cayley tree to be oriented towards the 
root node, as depicted in Figure [2] on page[6j thus, for an edge (w, v), v is closer 
to the root node and therefore the parent of w. The number of Cayley trees 
of size n is T n = n n ~ 1 . This formula is attributed to Arthur Cayley and thus 
also referred to as “Cayley’s formula”. The exponential generating function 
T( z ) := E n >o z of Cayley trees, the so-called tree function, is defined by 
the functional equation T(z) = ze T ^ z \ The tree function is closely related to the 
so-called (Lambert) W-function [5j that is defined by the functional equation 
z = W(z)e w ( z \ The simple connection between the Cayley tree function and 
the Lambert W-function is: T(z) = —W(—z). 

Definition 2. A function f from the set [n] into itself is called an n-mapping. 
Given an n-mapping f, its functional digraph Gf is defined to be the directed 
graph on the vertex set [n] where a directed edge ( i,j ) is drawn whenever f(i) = 

3- 

For an example of the functional digraph of a 19-mapping, see Figure [T| 
The structure of mappings is simple and is well described in HJ: the weakly 
connected components of their digraphs are simply cycles of Cayley trees. That 
is, each connected component consists of rooted labelled trees whose root nodes 
are connected by directed edges such that they form a cycle. In this context, 
we will call a node j that lies on a cycle in a mapping /, i.e., for which there 
exists a k > 1 such that f k (j) = j, a cyclic node. 

Using the symbolic method (see El for an introduction), this structural 
connection between Cayley trees and mappings can also easily be taken to the 
level of generating functions. Indeed, if T, C and M denote the combinatorial 
classes of Cayley trees, connected mappings and mappings and if T(z), C(z) 
and M(z) denote their exponential generating functions we have the following 
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connection: 

C = Cyc (T) =► C{z) = log ) , 

M = Set (C) =► M{z) = exp [C{zj) = • 

However, for the problem considered here the above relation between T and C 
cannot be applied directly. Instead we will use a decomposition of the objects 
with respect to the node with smallest label. 

Note that whenever we speak about a random n-mapping /, this means 
that / has been chosen uniformly at random among all n-mappings. Similarly, 
a random Cayley tree of size n is chosen according to the uniform distribution 
on all Cayley trees of the same size. 

In this paper, we shall generalize the concept of ascending runs from per¬ 
mutations to Cayley trees and to mappings. Let us start by defining ascents for 
permutations: An ascent of a permutation n is a position i for which it holds 
that 7 r(«) < n(i + 1). Similarly, a descent of a permutation n is a position i for 
which it holds that tv(i) > 7r (i + 1). The descents of a permutation 7r partition 
it into so-called ascending runs. These are increasing contiguous subsequences 
of maximal length. The following holds: if tv has k descents, it is the union of 
k + 1 ascending runs. 

Throughout this paper, we will write n— for the falling factorials n • (n — 
1 ) • • • (n — m + 1). 

3 A bijective proof of Cayley’s formula 

In [2], the present authors generalized the concept of parking functions to Cay¬ 
ley trees and to mappings. In this context, they presented a bijection between 
parking functions on trees and parking functions on mappings. The precise 
statement of this bijection can be found in Theorem 3.6 in [2]. When restrict¬ 
ing this bijection to trees and mappings only, we obtain (to the best of our 
knowledge) a new bijective proof of Cayley’s formula. We will be able to use 
this bijection in the following section that concerns ascending runs in trees and 
mappings. 

In the following, we will denote by T(v) the parent of node v in the tree T. 
That is, for v ^ root(T), T(v) is the unique node such that (■ v,T(v )) is an edge 
in T. For the sake of convenience, let us define T (root(T)) = root(T). 

Theorem 3 (Corollary of Theorem 3.6 in [2]). For each n > 1, there exists a 
bijection p from, the set of pairs ( T , w), with T £ T n a tree of size n and w £ T 
a node ofT, to the set of n-mappings. Thus 

n ■ T n = M n , for n > 1. 

Proof. Given a pair (T,w), we consider the unique path w root(T) from the 
node w to the root of T. It consists of the nodes V\ = w, V 2 = T(v i), ..., Vt + \ = 
T(vi),... ,v r = root (T) for some r > 1. We denote by I = (ii,...,i t ), with 
i\ < *2 < • • • < h f° r some t > 1, the indices of the right-to-left maxima in the 
sequence v±, V2, ■ ■ ■, v r , i.e., 

i £ I Vi > Vj, for all j > i. 
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The corresponding set of nodes in the path w root(T) will be denoted by 
Vi := {vi : i £ /}. Of course, if follows from the definition that the root node is 
always contained in Vj, i.e., v r £ Vj. 

We can now describe the function p by constructing an n-mapping /. The 
t right-to-left maxima in the sequence v\, V 2 , ■ ■ ■, v r will give rise to t connected 
components in the functional digraph Gf. Moreover, the nodes on the path 
w root(T) in T will correspond to the cyclic nodes in Gf. We describe / by 
defining f(v) for all v £ [n], where we distinguish whether v £ V/ or not. 

(a) Case v ^ Vr. We set f(v) := T(v). 

( b ) Case v £ Vj: We set /(wq) := v\ and f(vi e ) := T for l > 1. 

This means that the nodes on the path w -w root(T) in T form t cycles 

Ci := Oi, .. .,v h ), C t := ...,v r = v it ) in Gf. 

It is now easy to describe the inverse function Given a mapping /, we 
sort the connected components of Gf in decreasing order of their largest cyclic 
elements. That is, if Gy consists of t connected components and Cj denotes 
the largest cyclic element in the «-th component, we have c\ > C 2 > ■ ■ ■ > Ct- 
Then, for every 1 < i < t, we remove the edge (cj,di) where di = f(ci). Next 
we reattach the components to each other by establishing the edges {a,d i+ i) 
for every 1 < i < t — 1. This leads to the tree T. Note that the node Ct is 
attached nowhere since it constitutes the root of T. Setting w = di, we obtain 
the preimage ( T,w ) of /. □ 

Example 4. Taking the image of the pair (T, 1) where the tree T is depicted 
in Figure [2] leads to the mapping described in Figure[l] We consider the unique 
path from the node labelled 1 to the root of T. It consists of the following 
nodes: 1, 7,11,17,4,14 and 10. Within this sequence, the right-to-left maxima 
are 17,14 and 10 which are marked by gray nodes in the figure. When creating 
the image of (T, 1) under the map ip, the edges (17,4) and (14,10) are removed 
and the edges (17,1), (14,4) and (10,10) are created. H 

The bijection described above has the property of preserving ascents as well 
as elements that do not lie at the beginning of an ascending run: 

Proposition 5. Consider the pair ( T, w ) where T is a Cayley tree and w is 
a node in T. Let f be the image of (T, w) under the bijection ip described in 
Theorem [3| Then the following holds for every node v in T: 

1. T(v) > v f(v) > v, 

2. There exists a node x < v in T with T{x) = v 

<*=> There exists an element x < v in f with f(x) = v. 

Remark 6. These properties are fulfilled because of the following observation: 
When creating the image of some (T, w) under the map <p, the only edges that 
are removed are decreasing ones. The edges that are created instead in Gf are 
also decreasing ones. 

Proof of Proposition [5| 
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Figure 2: Taking the image of the pair (T, 1) where T is depicted above leads 
to the mapping described in Figure |T| The unique path from 1 to the root is 
marked by dashed edges. Right-to-left-maxima on this path are marked by gray 
nodes. 

1. As in the proof of Theorem [3] we distinguish whether the node v in T is 
contained in Vj or not: If v is not contained in Vj, f(v) = T(v) and thus 
the statement is clearly true. If v € Vi, v is a right-to-left maximum in 
the unique path from w to the root of T. It thus holds v = Vi e for some 
£ and clearly T(y) < v. We need to show that this implies that f(v) < v 
holds as well. By definition, we have /(vq) := V\ and f(vi e ) := T (vi e _ 1 ) 
for l > 1 Thus, f(vi t ) is either equal to v ig or f(v ie ) is not a right-to left 
maximum in which case it is clearly strictly smaller than Vi f . This proves 
the first statement of this proposition. 

2. (=>) : If the node v has a child x with x < v, then x cannot belong to the 
set of nodes Vj. Thus f[x) = T{x) = v and in the mapping / the node v 
has a preimage x with x < v. 

{<=) '■ Let x be an element with x < v and f{x) = v. Suppose x = v lf for 
some l. Then v = f(x) = T(vi e i ) < x, a contradiction. Thus x £ Vj and 
we have f(x) = T(x). Therefore x is a child of v in T that fulfils x < v. 

□ 

4 Main results 

We now consider the notion of ascending runs in a mapping. Generalizing the 
definition of ascending runs for permutations, an ascending run in a mapping 
/ is a maximal ascending sequence i < f(i) < f 2 {i) < ■ ■■ < f k {i), he., an 
ascending sequence which is not contained in a longer such sequence. 
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i is the start of a run 


Figure 3: The label pattern studied in this paper: Counting ascending runs is 
equivalent to counting the occurrences of the label pattern depicted above. 



Figure 4: A 19-mapping with 13 ascending runs. If an element is the first 
element in a run, its node is gray; otherwise it is white. Increasing edges, i.e., 
edges (i, j) with i < j, are marked by dashed lines. 


Given a specific mapping / and its functional digraph G/, the number of 
runs in / can be counted easily by considering the first element of a run: Indeed, 
an element j is the first element of a run iff each label i of a preimage of j is not 
smaller than j , i.e., (/(*) = j) => (i > j). We thus need to count nodes in the 
mapping that form the label pattern depicted in Figure [3] Note that elements 
with no preimages at all, i.e., terminal nodes in the graph, are of course always 
at the beginning of a run. 

Example 7. Consider the 19-mapping depicted in Figure [4] In its functional 
digraph, the nodes that correspond to first elements of an ascending run are 
marked by gray nodes. There are 13 such nodes and thus our mapping has 13 
ascending runs. Moreover, edges (*, j) where i < j are marked with dashed lines 
and edges with j > i are marked with full lines. An ascending run is thus a 
path of maximal length consisting of dashed edges only (or no edges at all). An 
element that is the start of an ascending run can be characterized as a node 
where no incoming edges are dashed. H 

Recall that for permutations the following simple connection between de¬ 
scents and ascending runs holds: if a permutation 7r has k descents it can be 
decomposed into (fc + 1) ascending runs. Moreover, by taking the complement 
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Figure 5: Decomposition of a tree with respect to its smallest node. 


7 r c of 7 r, i.e. the permutation mapping i to n — n(i) + 1, one obtains a permuta¬ 
tion with k ascents. Thus, studying the distribution of the following parameters 
in permutations is equivalent: ascents, descents, ascending runs and descending 
runs. This is clearly not the case for mappings: The presence of k descents 
does not yield a decomposition into (k + 1) ascending runs. For example, in 
the mapping depicted in Figure [4] there are 10 ascents and thus 9 descents but 
the number of ascending runs is equal to 13. Thus, the study of ascending 
runs has to be performed separately and requires, as we will see, more involved 
techniques than the one of ascents. Concerning the study of ascents, recall that 
the probability of i being an ascent in a random n-mapping is simply equal to 
(n — i)/n. See the work of Egecioglu and Remmel [7i as well as Clark [4] for a 
treatment of ascents in mappings and Cayley trees. 

In the following, R n will denote the random variable that counts the num¬ 
ber of ascending runs in a random n-mapping. Because of the combinatorial 
structure of mappings, we will first analyse the corresponding random variable 
for trees before we study R n . 

4.1 Ascending runs in Cayley trees 

To study R n we first consider the corresponding quantity in Cayley trees. Thus, 
we introduce the random variable F n which counts the number of runs in a 
random Cayley tree of size n. Let us further introduce the generating function 

z n z n 

n*,«) = EE w = = 

n> 1 m> 0 n> 1 m> 0 

where T n = n n_1 is the number of trees of size n and F n rn is the number of 
trees of size n with exactly m ascending runs. 

Clearly, m = S m ^. In order to establish a recurrence relation for the 
numbers F n rrl . we decompose a tree of size n > 2 with respect to the node 
labelled 1. Two different cases, as shown in Figure [5j have to be considered: 

Case (1): The node 1 is the root node. In this case a set of r > 1 subtrees 
A, T 2 ,..., T r is attached to the root 1. Since the root has the smallest of all 
labels, it will always constitute a run of its own. Thus, if the subtrees Ti have 
respective sizes rii and m t runs each, the entire tree will consist of n nodes and 
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have to runs iff ^[ =1 n* = n — 1 and m i = to — 1. The contribution to 

F nrn is therefore: 


pi 1 ) = 

n,m 


E 

r>l 


= £ 


F, 


n i ,mi 




yy m=n— 1, 

y: 



where the multinomial coefficient reflects the number of possibilities of redis¬ 
tributing the labels in [n — 1] order-isomorphically to the subtrees. 

Case (2): The node 1 is not the root node. In this case the node I is attached 
to some node labelled i of a tree To of size no and with too runs and has itself 
a set of r > 0 subtrees T\, T 2 ,..., T r attached to it. Depending on where in 
T 0 the node I is attached, a new run is created or not. If 1 is attached at the 
beginning of a run in To no new run is created - there are mo possibilities of 
attaching 1 in this way. If 1 is however attached somewhere in the middle of a 
run, a new run will be created - for this case there are (n o — too) possibilities. 
The contribution to F n ^ m is therefore: 


Fn,m — y", m o ■ y ' y ' F n0>mo ■ ■ ■ F nrm: 

l<no<n—l, r> 0 ^2m—n—n 0 — 1, 

l<mo<m yy rrii—m—mo 



for the case that the node 1 is attached to the beginning of a run in Tq and 


Fn,rn ~ ^ ^ ( n 0 m o) * ^ ^ ^ -^n 0 ,m 0 * * * ^n r ,m. 

l<no<n—1, r>0 yy m=n—no — l, 

l<m Q <m yy rrii=m—mo — l 



in the other cases. 

The goal is now to translate this recurrence relation for T„ jm = Fnjn + 
F^m + Fi 2 m into the language of generating functions. In order to alleviate 
notation we write F for F(z, v), F z and F v for the partial derivatives of F with 
respect to z and v. 

First, let us note that for Fn]m we have the following: 

F ^ = E^7T^ E [z ni v^]F---[z n ^]F 

r> 1 yy m=n— 1, 

yy rrii—m— 1 

= (n — l)\[z n ~ 1 v m ~ 1 ] exp(T) = (n - l)\[z n v m }zvey^>{F). 

Recall that to 0 • T no>mo /(n 0 )! = [z no v m °]vF v . We then obtain for Fn,m- 

F%$ = (n- 1)! E {[ zn °v rn °}vF v ) • ([^ n_ " 0_1 n m_m °] exp(T’)) 
n 0 ,m 0 

= (n - 1 )\[z n v m ]zvF v exp(T). 

Similarly, for F^m'. 

Fg = (n - 1 )\[z n v m }zv ( zF z - vF v ) exp(F). 
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Thus, we can write F n rn as follows: 

F n , m = n\[z n v m ]F = (n - l)\[z n v m ]zv exp(F) (1 + + (1 - v)F v ), 

or equivalently, 

EE n[z n v m }F = zF z = zwexp(F) (1 + 2 .F 2 + (1 - v)F v ). 

n> 1 m>l 


Thus the bivariate generating function F(z, v) is defined by the following first- 
order quasi-linear partial differential equation: 

(l - zve F{z ' v F z (z,v) = v(l - v)e F(z ’ v) F v {z,v) + ve F{z ' v \ 


with the initial condition F(0,v) = 0. 

This PDE can be solved using the method of characteristics (see, e.g., ID- 
Let us regard z, v and F as variables depending on t. We then obtain the 
following system of ordinary differential equations, the so-called system of char¬ 
acteristic differential equations: 


dz 

dt 

dv 

dt 

dF 

dt 


= 1 -zve F ( z ' v \ 

= e F{z ’ v) -v(v- 1), 
= ve F ^ v \ 


(la) 

(lb) 

(lc) 


We are now looking for first integrals of the system ([I]), i.e., for functions 
C(z, v, F) that are constant along any solution curve - a so-called characteristic 
curve - of 0 - From (|lb[) and (flc|) one obtains the following differential equation 


dv 

dF= V 1 


which has the general solution 


v = Ci • e F + 1. 


Thus a first integral of ([T]) is given by: 

Ci (z,v,F) =a = 


„-F 


(v ~ !)• 


( 2 ) 


From (la) and (lc) one obtains, after substituting v = Ci ■ e F + 1, the following 
differential equation 

dz 1 

dF (cie F + 1) • e F 

The general solution of this first order linear differential equation is 

F — In (ci • e F + l) + c 2 


After backsubstituting ci = e F • (v — 1) one obtains the following first integral 
of |l]) which is independent of 0: 

C 2 (z, v, F) = c 2 = z ■ e F - F + ln(v). (3) 
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The general solution of ([I]) is thus given by: 

G(Ci(z,v,F),^ 2 (z,v,F)) = const., 


where G is an arbitrary differentiable function in two variables. One can also 
solve this equation with respect to the variable z and obtains: 

z = e~ F ■ ^F-ln(u) +g ^ , 


where g is an arbitrary differentiable function in one variable. In order to specify 
this function g , we plug the initial condition F( 0, v) = 0 into this equation. This 
leads to: 

g (v - 1) = ln(u), 

and finally we obtain that the solution of ([I]) is given by the following functional 
equation: 


z = 


In ( e F(, '"*- 1+ » ) 


e F(z,v) 


(4) 


4.2 Ascending runs in (connected) mappings 

Having computed the generating function F(z,v ) of ascending runs in trees in 
the previous section, we now turn to our actual problem, namely the enumer¬ 
ation of mappings by size and number of ascending runs. We introduce the 
following bivariate generating function 

z n 

R (z,v) = v m , (5) 

n>0 m >0 

with M n = n n the number of n-mappings. Actually, it is easier to study the 
quantity for connected mappings first. We thus introduce the corresponding 
generating function C(z, v). Due to the combinatorial construction of mappings 
as sets of connected mappings, it holds that 

R(z,v) =e G( - z ' v \ 

Note however that such a simple relation does not hold when passing from 
Cayley trees to cycles of Cayley trees. Indeed, the number of ascending runs 
in a connected mapping is not equal to the sum of the numbers of ascending 
runs in the Cayley trees constituting this connected mapping. Thus, it does not 
hold that C(z,'v) is equal to log (1/(1 — F(z,v))) as one might be tempted to 
think. In order to establish a functional equation for C(z,v), we will therefore 
decompose a connected mapping graph according to the node labelled 1 as done 
before for Cayley trees. Here, we have to consider three different cases, two of 
which are depicted in Figure [6] 

Case (1): The node 1 is a cyclic node in a cycle of length one. In this case the 
mapping graph is simply a tree with root node 1 and an additional loop-edge 
(1,1). We thus have exactly the same situation as in the first case for trees. See 
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Figure 6: Decomposition of a connected mapping with respect to its smallest 
node. 


the left hand side of Figure [5j The contribution to C n , m is therefore: 



E 


F n 


r>l y m=n— 1, 
y mi=m— 1 

= (n-l)\[z n v m }zvexp(F). 



Case (2): The node 1 is a cyclic node in a cycle containing more than one 
element; this case is depicted on the left hand side of Figure[6] The structure of 
such a mapping can be understood as follows: There is a (non-empty) tree To 
in which the node 1 has been attached to some node i. Moreover, the root of To 
is attached to the node 1, thus forming a cycle consisting of the path from i to 
the root of To and i. Since 1 can have arbitrarily many children, we pick a set of 
r > 0 subtrees and attach them to 1. Whether the node 1 contributes to a new 
run or not depends on whether it is attached at the beginning or somewhere in 
the middle of a run (as in the second case for trees). Indeed, if To is of size no 
and has mo runs itself, there are mo possibilities of attaching 1 at the beginning 
of a run in To thus not creating a new run. Moreover, there are (no — m o) 
possibilities of attaching 1 somewhere in the middle of a run and of creating a 
new run. The contribution to C n , m is therefore: 


C n,m — 'y ' m 0 y ' j y ] Fn 0 ,m 0 ' ' ' Fn r ,m r 

l<n 0 <n—l, r >0 JZ rii—n— n 0 — 1, 
l<m 0 <m y mj—m—mo 



H - 5> m 0 ) ' ^ ^ j ^ ^ ^no,mo * * * Fn r ,m r 

l<no<n— 1, r>0 ^nj=n-no — 1, 

l<mo<m y rrii—m—mQ — l 



(n - l)![z n ^ m ]^exp(F) (zF z + (1 - v)F v ). 


Case (3): The node 1 is not a cyclic node; this case is depicted on the right hand 
side of Figure [6j Here, the node 1 is attached to some node i of a connected 
mapping Co and a set of r > 0 trees are attached to 1. Again, the node 1 
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contributes to a new run or not depending on where it is attached. Similarly as 
in the second case, we obtain that the contribution to C n , m is: 


c( 3 ) — \ ( ^ no ’ mo p ... p 

'~'n,m — / y '"-0 / y | / y r ni,mi r n r ; 

l< n 0 < n— r“>0 22 ni—n—no — 1, 

l<m 0 <m '22 mi=m—mo 



+ J^( U ° ~ W °) ' C T° H F ni>Tni ■ ■ ■ F„ r> . 

l<n 0 <n—1, r>0 22 rii—n—nQ — l. 

l<mo<m mj—m—mo — 1 



= (n - l)![z n w m ] 2 t>exp(F) ( 2 C 2 + (1 - v)C„). 

Adding up the contributions of these three cases and summing over n > 1 
and 77i > 1 we obtain the following: 

C z = vexp(F) • (1 + (1 - v)F v + zF z + (1 - v)C v + zC z ). 

Thus, the bivariate generating function C(z,v) is defined by the following first- 
order linear partial differential equation: 


C z (1 zv • exp(F)) - v(l - v ) exp (F)C V = F z , 

where F(z,v) is defined by equation Q. 

Using the auxiliary function F[(z, v) that is defined by 

(H) = ^PiF)-l+v 

V 

and that satisfies the following functional equation 

H 

Z Tt ' Z ; 

ve n + 1 — v 

we can solve this PDE. We find that the solution is given as follows: 


C(z, v) = In 


0 H(z 


+ 1 — 7 ) 


ve H ( z ’ v ) (l — H{z , w)) + 1 — ■ 


( 6 ) 

(7) 


with H(z, v) as given above. Checking that this function indeed is a solution to 
the partial differential equation can be done easily by using implicit differenti¬ 
ation in order to compute the partial derivatives. 

We thus also obtain the solution for ascending runs in arbitrary mappings: 


R(z, v) = e 


- p C(z,v) _ 


e H{z,v) + ! _ v 


ve H(z,v) — H(z,v)) + 1 — V 

H/z 1 


H/z — Hvexp(H) \ — zv ■ ex^F[{z,v))' 


with H(z,v) satisfying 0- 


( 8 ) 
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4.3 Exact enumeration formulae 

Having found the generating functions F(z,v) and R(z,v) for trees and map¬ 
pings, respectively, we can prove the following interesting connection between 
ascending runs in trees and in mappings: 

Theorem 8. Let F nm denote the number of Cayley trees of size n with exactly 
m ascending runs and R n ,m the number of n-mappings with exactly m ascending 
runs. Then for all n > 1 and m > 1 the following identity holds: 

Rn,m ^ ' Rn,m- 


Proof. We will prove the assertion above at the level of generating functions. 
Combining equation § and ([8]) we have: 

R(z,v) = — 


1 — z - (exp(F(z, v) — 1 + v)) 
Moreover, implicit differentiation of Q with respect to 2 leads to: 

= exp(F(z,*;))-l + i> 

1 — 2 • (exp(F( 2 , v) — 1 + v)) 

Thus one easily sees that 

R(z,v ) = 1 + 2 • .FI ( 2 , v). 

Extracting coefficients for n > 1 and m > 1 leads to 

Rn.m = n\[z n v m ]R(z,v) = n\[z n ~ 1 v m ]F z (z,v) = n ■ 


□ 


A combinatorial proof of this result can be found in the following section. 
Combining this result with the functional equation for the generating func¬ 
tion F(z,v) we obtain the following enumeration formulae: 

Theorem 9. The number F n ^ m of Cayley trees of size n with exactly m runs 
is given as: 

F n , m = (n- 1)^4 

[rri J 

and the number R n ,m of n-mappings with exactly m runs is given as: 



(9) 


where {^} denotes the Stirling numbers of the second kind that count partitions 
of the set [n] into m non-empty subsets. 
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Proof. Extracting coefficients from the functional equation ([4]) for F(z,v), the 
generating function of ascending runs in Cayley trees, can be done by applying 
Lagrange’s inversion formula (see, e.g., [TT] ) with tp(F) = exp(E)F/ln 
and twice Cauchy’s integral formula: 




1 , e nF • F r ‘ 

x r —n-li c 1 


( \ ? 

e ~^) 


l 


n ■ 27 xi 
1 

n ■ 2iri 


0 nF 




rdF 


In 


(■ e s ■ v — 1 + v) 


n 1 • e s v 


dS 


= ■ v - 1 + v) 

n— 1 


n ~ 1 ■e s v 


= - y_1_y 

n ^ (n — 1 — fc)! 7 —^ 
fc= o v ' t=o 


. - 1\ tv £ 


k\ 


(v - 1) 


i-i-e 


= -y 

n ' 


£=0 


n — 1 


v e (v — 1) 


l-l -l (l + 1) 


n— 1 


(n- 1)! 


When passing from the second to the third line above we use the substitution 
S = In 


(^) 


for which it holds that dR = - — -k^-dS = . f e , , - dS. 


ve s -l-\-v 


What remains to do now is to extract the coefficients in v: 

n-1 , 

F n}m =n\[v m z n }F(z,v) = Yl { e J^+ir-'y-^Kl-y 


n —1 

= £ 

£=0 


n — 1 




m—i— 1 


Finally, we obtain the following explicit solution for the numbers F n>m = 
T n ¥{F n = m} = n\[z n v m ]F(z, v) of size-n trees with exactly to runs: 

v 7 e=o v 7 

These numbers can also be described with the help of the Stirling numbers of 
the second kind. Since the Stirling numbers of second kind satisfy 


TO TO! 


e-o 




( 10 ) 


we obtain the following for n > 1: 
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= («- !)—{”}• (H) 

Using the connection between the numbers F nrn and R n ,m as obtained in The¬ 
orem [HJ we obtain the corresponding result for runs in mappings. □ 

4.4 Bijective proofs 

In the previous section, we proved two results about ascending runs in trees and 
mappings where a combinatorial, i.e., bijective explanation is desirable. First, 
the fact that R n ^ m = n ■ F nm and second, the fact that the numbers R n , m can 
very easily be expressed with the help of the Stirling numbers of the second 
kind. In this section we will thus provide bijective proofs for these statements. 

Theorem 10. There is a bijection between pairs (T,w), where T is a tree of 
size n with exactly m ascending runs and w is a node in T and n-mappings f 
with exactly m ascending runs. 

Proof. This statement follows quite directly from Proposition [5] For a pair 
( T,w ), let / be its image under the bijection described in Theorem [3] We need 
to show the following: A node v is at the start of a run in T if and only if v 
is at the start of a run in /. As remarked earlier, a node v is at the start of 
a run in a tree or mapping if and only if all children or preimages x of v have 
labels larger than x. Then the statement of this theorem follows by negating 
the second property in Proposition [5j □ 

Example 11. As described in Section [3j the pair (T, 1) where T is depicted in 
Figure [2] is identified with the mapping / depicted in Figure [IJ As can be seen 
very easily both T and / have 13 ascending runs and the nodes that are at the 
start of a run are the following: 1, 2, 3,4,5,6, 8, 9,12,13,15,18 and 19. H 

Now let us turn to the second interesting fact that asks for a bijective proof, 
namely that the number R n ^ m of n-mappings with m ascending runs can be ex¬ 
pressed with the help of the Stirling numbers of the second kind as in Theorem[9] 
We shall prove the following: 

Theorem 12. There is a bijection between the set of n-mappings with exactly 
m runs and the set of pairs (S,x), where S is a set-partition of [n] into m parts 
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Cq 

II 

^0 

ni = 13 

S 2 = {18} 

n 2 = 13 

S 3 = {17,13} 

n 3 = 1 

S 4 = {16,9} 

714 = 4 

S 5 = {15} 

n 5 = 7 

S 6 = {14,4} 

n 6 = 4 

S 7 = { 12 } 

5 

-vl 

II 

-a 

S 8 = { 11,7,6} 

n 8 = 17 

S 9 — { 10 , 8 } 

ng = 10 

S 10 = {5} 

n w = 17 

II 

to 

nn = 10 

S 12 = { 2 } 

m 2 = 1 

Si 3 = { 1 } 

ni3 = 7 



Figure 7: Example of the bijection described in the proof of Theorem [12] for the 
mapping depicted in Figure [T] 


and x = (ni ,... ,n m ) is an integer sequence of length m. The set partition is 
given as S = (Si, S 2 , ■ ■ ■, S m ) where the parts are ordered decreasingly according 
to the largest element in each part, i.e., it holds max(5i) > max(52) > ••• > 
rnax(«S' m ). The sequence x then has to fulfil the following restriction: Uj £ 

N \ (Ui=i min{£ £ Si : i > max(S'j)}). 

Remark 13. The idea of the bijection is to successively decompose the mapping 
into ascending runs. This is done by starting with a run ending at the largest 
element of the mapping, then one ending at the next-largest element that has 
not been involved yet, and so on. The runs then correspond to blocks of the 
partition. In order to keep track of how these runs where “glued” together and 
to be able to reconstruct the mapping, we additionally store the image of the 
last element of each run in the sequence x. 

Proof of Theorem\l First we remark that this indeed will prove that 



since the number of set-partitions of [n] into m parts is given by {^} and the 
number of sequences x satisfying the restrictions is given by n ■ (n — 1) • • • (n — 
m + 1) = n—. 

To prove the theorem we consider an n-mapping with exactly m runs and 
iterate the following procedure, where we colour the elements of the mapping 
until all elements are coloured. 

• In the j-th step we consider the largest element in the mapping which has 

(7) 

not been coloured so far; let us denote it by s{ . Consider all preimages 
of 4 f) with a label smaller than and, if there are such ones, take the 
one with largest label amongst them; let us denote this element by sf . 
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(i) (?) 

Then iterate this step with si,', i.e., amongst all preimages of sf with a 

label smaller than sf take the one with largest label, which is denoted 

by After finitely many steps we arrive at an element , which 

does not have preimages with a smaller label. We then define the set 

Sj : = { S ^K ■ ■ ■ > s k )}• Note that in the mapping graph this corresponds to 

a path s]f. —)••••—> s^ —>■ s{' with increasing labels on it. 

• Additionally we store in rij the image of s^. Clearly s^ is in [n]. Due 
to the construction further restrictions hold: Indeed, if i < j, nj cannot 
be the smallest element in Si larger than s^ \ since otherwise s^ would 
have been chosen during the construction of the set Si. 

• Finally colour all elements of the mapping contained in Sj. 

Since the mapping contains exactly m runs and the smallest element in each 
set Sj corresponds to the minimal element of a run the procedure stops after 
exactly m steps. It thus defines a pair of a set partition S = (Si, ..., S m ) and 
a sequence x = (n±, ..., n m ) with the given restrictions. 

If the pair (S, x) is given, the corresponding mapping can easily be recon¬ 
structed. Indeed, the partition S gives us a decomposition of the mapping into 
ascending runs and the sequence x tells us how these runs have to be linked to 
each other. The inverse of this bijection can therefore be defined in a straight¬ 
forward way. □ 

Example 14. The construction of the partition S and the sequence x for the 
mapping described in Figure [l] can be found in Figure [7J Let us exemplarily 
explain how the set S§ is constructed. At this point, the elements in 
i.e., 19,18,17,13,16,9,15,14,4 and 12, have already been coloured. Thus, the 
largest element that has not been coloured so far is 11 = s{ ; . For S 2 , we 
consider the preimages of 11 that have a label smaller than 11. The only such 
element is 7 and thus s® = 7. Next, the preimages of 7 are 6 ,12 and 15 and thus 
S 3 = 6 . Since 6 does not have any preimages, we stop here and S$ = {11, 7, 6 }. 
Since the image of 11 is 17, we set rig = 17. H 

4.5 Limiting distribution 

In order to show that the random variable R n has linear mean and variance and 
that it is asymptotically normally distributed when suitably standardized, we 
apply H.-K. Hwang’s quasi-power theorem: 

Theorem 15 (Hwang’s quasi-power theorem |T3]). Let (A n )„ e pj be a sequence 
of discrete random variables supported by No- Assume that the moment gener¬ 
ating functions g n (s ) = E ( e sXn ) = = Li)e sk are analytic in a fixed 

complex neighbourhood of s = 0 in which they have an expansion of the form 

9 n(s) = exp {<t> n U{s) + V{s)) • ^1 + 0 > 

where 4> n ,K n —> 00 and U(s), E(s) do not depend of n and are analytic at 
s = 0. Assume finally that U(s) satisfies the so-called “variability condition” 
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U"( 0) ^ 0. Under these conditions, the mean and variance of X n satisfy 

E(X n ) = <j> n U'( 0) + V'(0) + O (O. 

V(x n ) = ^c/"(0) + y"(0) + 0(0- 


Furthermore, X n is after standardization asymptotically normally distributed, 
with speed of convergence 0(K~ k + (fn 1 ^ 2 ): 


X n - E(X„) \ 

y/WQ 


<f>(x) + 0 




where <f>(:r) is the distribution function of the standard normal distribution. 
We then obtain the following: 


Theorem 16. Let R n denote the random variable counting the number of as¬ 
cending runs in a random n-mapping. Then the expectation and variance satisfy 

E (R n ) = (1 — e -1 )n + 0(1) and 
Y(R n ) = (e- 1 - 2e~ 2 )n + 0(1). 


Furthermore, the standardized random variable 

R n - E (Rn) 

VWQ 

is asymptotically Gaussian distributed with a rate of convergence of order 0(n~ x ! 2 ). 


Remark 17. Note that Theorem [8] implies that the random variable R n count¬ 
ing the number of runs in ro-mappings and the random variable F n counting 
the number of runs in trees of size n follow exactly the same distribution. The 
theorem above thus holds for F n as well. 

Proof of Theorem \lb\ Before we get our hands on the bivariate generating func¬ 
tion R(z,v ) defined in equation (pH), we will need to take a closer look at the 
function F[(z,v) from equation that is involved in the functional equa¬ 
tion ([8]) for R(z,v). Once we have a singular expansion of H(z,v), we can 
obtain an expansion of R(z,v). Extracting coefficients [z n ] in R(z,v ) will 
then allow us to obtain an asymptotic expansion of the probability generat¬ 
ing function p n (v) = = fc)u fe and the moment generating function 

9n(s) = p n (e s ). This will finally allow us to apply the quasi-power theorem and 
to obtain the statement of this Theorem. 

Asymptotic expansion of Ft. The function Ft = H(z,v ) is defined by a func¬ 
tional equation of the form Ft = zip(H) where the function ip is defined as 
p(u) = v ■ e u + 1 — v. Singularities of Ft can be found with the help of singu¬ 
lar inversion (see Theorem VI.6 in USD- For this, we need the characteristic 
equation of ip: 

tp(r) — r • p'(r) = 0 
<t=> ve T + 1 — v — t ■ ve T = 0 

1 — v , , 

•£=> r = 1 +-. (12) 

ve T 
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For v = 1 we have r = 1 and whenever v is close to 1 a unique solution to 
equation (12) exists as can be seen by applying the Implicit Function Theorem. 
Note that r is a function of v. This will however be omitted in the notation and 
in the following, r will always denote the unique solution of equation (12). The 
unique dominant singularity of H 1 considered as a function in z, is at z = p , 
with: 


1 


1 


P = 


<p'M 


d 


1 —V ’ 


if v = 1 , 
if v 7 ^ 1 . 


(13) 

(14) 


Before we continue by giving the asymptotic expansion of iJ, we want to charac¬ 
terize its unique dominant singularity p with the help of a functional equation. 
From (12) and p = t/<p(t) it follows that: 

t = 1 + p{ 1 — v) and 

r = p ■ (ye T + 1 — v) = p ■ ^ue 1+p(1_ ^ + 1 — v'j . 

Dividing by p then leads to 
1 


+ 1 — v = ve 


- „D+p( i-D 


1 - v, 


or, equivalently, 


Thus 


1 


P = 


e ■ ePD-D' 


(1 — v)pe 


Pi i-v) = 1 ~ v 
v • e 


(15) 


and p can be expressed with the help of the Lambert VF-function for v ^ 1 but 
in a neighbourhood of 1 : 


1 


P = 


1 — v 


-W 


1 — v 


Note that the above expression is not defined for u = 1. However, taking 
limits as v tends to 1 , the right-hand side tends to exp(— 1 ) as expected from 
equation (14). 

The function H admits the following asymptotic expansion around its unique 
dominant singularity p, i.e, for z —> p: 


H{z,v) = t — 


' Mr) 

<P"(t) 


l--+K-[l--)+o{l-- p j 

3/2 


A 3/2 




where K is a computable constant that will not be needed explicitly in the 
following. 
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Asymptotic expansion of R. We can now use this result in order to describe 
the analytic behaviour of the function R(z,v). Indeed, R inherits a unique 
dominant singularity at z = p from H. This can be seen as follows: Let us 
first recall that (R(z,v))~ 1 = 1 — zv ■ exp (H(z,v)). Thus, R has a singularity 
whenever zv • exp (H) = 1 or equivalently, whenever H = — ln(zn). Using the 
functional equation for H given in equation 0. this is equivalent to: 


— ln(n>) 

v ■ exp(— ln(zu)) + 1 — v 


& 1 


—ln(zn) 

1 + 2 — zv 


2 exp ( 2(1 


v)) 


1 

ev' 


which is exactly the same as equation (151 for p. Thus z = p is the unique 
dominant singularity of R. 

In order to provide the asymptotic expansion of R for 2 —> p, let us write 2 
as (2 — p) + p and H as (H — r) + r: 


t = 1 - vpe T ■ e H ~ T + vp ( 1 - ^ ) e T ■ e n ~ T = 1 - e n ~ T + ( 1 - y 


H-i 


H-7 


H-7 


since vpe r = 1 by definition of p (see equation @). Using the power series 
expansion of the exponential function and the asymptotic expansion for H(z , v) 
we obtain the following (for 2 —> p): 


„H-- 


and thus 


= 1 + (H-t) + 0((H - t) 2 ) = 1 - VzP^Jl- - p + o(l - 

i = ^yUf +o(l _£). 


Finally, with r = 1 + p(l — v), we obtain 
R(z,v) = 


1 

y/2y/l+p(l-v)Jl^j 


1 + 0 [J1-- 


( 16 ) 


Asymptotic expansion of p n (v). By a standard application of singularity 
analysis (see m for an introduction) we obtain asymptotics of the coefficients 
of R(z, v ): 


[z n ]R(z,v) 


[z n ]R(zp, v) 

1 

V2\/l + p{ 1 ^ v) 

1 

+ pC 1 - v ) 
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Since p n (v) = n\/n n • [z n ]R(z, v), we obtain the following expansion of the 
probability generating function by applying Stirling’s formula: 


Pn{v) = 


1 


1 


n\ 1 

n n p n y/2yjl + p( 1 - v) V™ 
1 \\ 1 


1+0 


1 + o 

e' u n u \ \n J J p" 

1 = --L- (l +o(~^= 

+ p{l - v) V™ V W™ 


1 


1 


(e • p) n ■ sjl + p(l - v ) 


1 + 0 


1 


(17) 


Applying the quasi-power theorem. Expansion ( |17[ ) now immediately gives 
an asymptotic expansion for the moment generating function: 


9n(s) = 


1 


1 + 0 ( 4 = 

n 


(e • p) n ■ yjl + p(l - e s ) 

= exp (-71 • ln(e ■ p) - In (y/l + p(l - e s )^ • ^1 + 0 . 

This is precisely the situation where the quasi-power theorem, Theorem |15[ can 
be applied. The involved functions are here defined as follows: 


Kn = \fn 


U(s) = — ln(e • p), 

and V(s) = — In (^/l + p( 1 — e s )^j . 


Note once again that p is not a constant but depends on e s . We only need 
to check the variability condition U"( 0) ^ 0; the other conditions are clearly 


fulfilled. For this purpose, let us use equation (13) and write U(s) as follows: 


C/(s) = ' lll (^) = - 1 + s + T - 

We thus obtain that the first and second derivative of U(s ) are given as follows: 

U'(s) = 1 + T\e s )-e s , 

U"{s) = e 8 ■ (r"(e s ) + t' (e s )), 

where r is considered as a function in v = e s . In order to determine t' and t” . 
we will make use of implicit differentiation and the functional equation for r 


given in (12): 


1 — v 


t\ implying t' = 


1 


Furthermore, 


v ■ (v — 1 — ve T ) 

(v — 1 — ve T ) + v ■ (1 — e T — ve T ■ t') 
v 2 (v — 1 — ve T ) 2 


+ 


2ve T + 1 — 2v 


v(v — 1 — t>e T ) 3 v 2 (v — 1 — ve T ) 
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Recalling that r = 1 for e s = v = 1, we obtain: 

C/'(0) = 1 + t'(1) = 1 + = 1 - e" 1 w 0.632..., 

—e T '- L > 

U"( 0) = r"( 1) + r'(l) = + 2 ^ L - e _1 = -2e“ 2 + e” 1 « 0.0972 .... 

—e 6 e z 

In a similar way, we could also determine V^O) and V"(0). Since the calcu¬ 
lations are rather tedious and the additive constants in E(i?„) and Y(R n ) are 
not of particular interest, we omit this here. It suffices to note that F'(0) and 
V"(Q) are both real constants. This finally gives us the desired result that 

E(i? n ) = (1 — e -1 )n + 0(1) and Y(R n ) = (e -1 — 2 e~ 2 )n + 0(1), 

and that the standardized random variable ( R n — E (R n ))/\/Y(Rn) converges 
to a standard normal distribution. Moreover, since k” 1 = n^ 1 / 2 = 4> n 1 ^' 2 , the 
speed of convergence is of order 0(?r -1 / 2 ). □ 

5 Summary of the results and future work 

Let us briefly discuss the results of this paper. 

For the number of ascending runs in Cayley trees and mappings, we obtained 
exact enumeration formulae and a connection to the Stirling numbers of the 
second kind that we could also explain with the help of bijective proofs. We 
also obtained limiting distribution results and showed convergence to a normal 
distribution. Let us compare this to the known results for the Eulerian numbers 
that enumerate permutations by their length and their number of ascents. 

The random variable X n counting the number of ascents in a random permu¬ 
tation of length n satisfies a central limit theorem (3l. Indeed, with E(A'„) = 
and Y(X n ) = it holds that the standardized random variable (X n — 
K(X n ))/y/Y(X n ) converges in distribution to a standard normal distribution. 
Equivalently, the random variable Y n = X n + 1 counting the number of ascend¬ 
ing runs in a random permutation of length n satisfies a central limit theorem 
with E(y„) = and Y(Y n ) = Y(X n ). Thus, both for ascents and for ascend¬ 
ing runs, it holds that E ~ n/2 and Y ~ n/ 12 = 0.083 • n. Here, we could show 
similar results for ascending runs in trees and mappings. Indeed, we also ob¬ 
tained a linear mean and variance and convergence in distribution to a standard 
normal distribution of the standardized random variable. However, the involved 
coefficients are not the same: 

E(i? n ) ~ (1 — e _1 )n with 1 — W 1 ss 0.632... and 
Y (Rn) ~ (W 1 — 2 e~ 2 )n with W 1 — 2e -2 « 0.0972_ 

The combinatorial decomposition of Cayley trees and mappings employed in 
this paper can also be used to study other problems in connection with label pat¬ 
terns. First, we remark that with the generating functions approach it is indeed 
possible to carry out a joint study of ascending and descending runs, although 
(as might be expected) the occurring functions are more involved. Furthermore, 
our method can also be applied to count the number of local minima or maxima, 
thus generalizing the notion of valleys and peaks in permutations. However, the 
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results contained in Okoth’s thesis [T5] are more general and thus the study of 
local minima and maxima was not included here. Moreover, preliminary stud¬ 
ies show that our method can also be applied to what we call up-ups. For an 
n- mapping /, these are elements i in [n] for which it holds that i < f(i) < f 2 (i). 
The notion of up-ups thus generalizes consecutive 123-patterns in permutations. 
Permutations avoiding this pattern have been studied, e.g., in [Bj. Finally, we 
remark that with the same techniques one could also attack corresponding label 
pattern problems for other combinatorial tree families, such as labelled ordered 
trees and labelled binary trees. 

All label patterns mentioned so far in this paper were consecutive patterns, 
i.e., they concern the images or preimages of nodes in a mapping. A more 
involved task would be to study non-consecutive label patterns in mappings. For 
instance, the simplest non-trivial non-consecutive pattern 21 corresponds to an 
inversion in the mapping. The number of inversions in trees has been studied in 
several papers, see, e.g., H2E3, and it seems that limiting distribution results 
carry over from Cayley trees to mappings. Of course, it would be interesting to 
obtain results also for longer non-consecutive patterns in mappings or Cayley 
trees. 
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